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ABSTRACT 

We seek to clarify the nature of running penumbral (RP) waves: are they chro- 
mospheric trans-sunspot waves or a visual pattern of upward-propagating waves? 
Full Stokes spectropolarimetric time series of the photospheric Si 1 10827 A line 
and the chromospheric He 1 10830 A multiplet were inverted using a Milne- 
Eddington atmosphere. Spatial pixels were paired between the outer um- 
bral/inner penumbral photosphere and the penumbral chromosphere using incli- 
nations retrieved by the inversion and the dual-height pairings of line-of-sight ve- 
locity time series were studied for signatures of wave propagation using a Fourier 
phase difference analysis. The dispersion relation for radiatively cooling acoustic 
waves, modified to incorporate an inclined propagation direction, fits well the 
observed phase differences between the pairs of photospheric and chromospheric 
pixels. We have thus demonstrated that RP waves are in effect low-/? slow-mode 
waves propagating along the magnetic field. 

Subject headings: Sun: infrared - Sun: magnetic fields - Sun: sunspots - Tech- 
niques: polarimetric - Waves 



1. INTRODUCTION 



The term running penumbral (RP) wave was created by IZirin fc Steinl ( 119721 ) to de- 
scribe chromospheric Ha velocity and intensity fronts that were observed moving out through 
sunspot penumbrae. Since then, a hos t of work has been carried out on re porting their prop 



erties (see, e.g., the series of papers by lChristopoulou et al.ll2000i . 1200 ll and iGeorgakilas et al. 
2000l ). while their exact nature has remained unidentified. Currently, the two most likely 
possibilities for the form of these oscillatory disturbances are: 



1. trans-sunspot waves generated in the umbra (e.g., by umbral flashes) and limited to 
the chromospheric layer, 
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2. a "visual pattern" resulting from field-aligned waves propagating up from the photo- 
sphere. 



To date, many findings point toward RP waves being due to the " visual pattern" scenari o 
(for an extensive discussion of this t opic see the recent r eview by iBogdan fc Judge! 120061 ). 
This seems especially likely now that ICenteno et al.l (120061 ) have successfully identified chro- 
mospheric 3-mi n umbral oscillat i ons a s prop agating, field-aligned, acoustic waves. However, 
recent work by iTziotziou et al.l (120061 . 120071 ) has not been able to decide between either of 
the two possible RP wave scenarios. Thus, the sum of the evidence is still not conclusive. 

In this paper we use velocity time series observations that possess a two-fold advantage 
over those of previous studies. The first is the simultaneous recording of photospheric and 
chromospheric lines, allowing the connection between velocities in the lower and upper at- 
mosphere to be accurately investigated. The second is the retrieval of the magnetic vector 
at both heights in the atmosphere through the use of full-Stokes spectropolarimetry, thus 
circumventing the need for any assumptions about the possible orientation of the magnetic 
field. 



OBSERVATIONAL DATA 



The dataset used here was obtained from the main sp ot of active region NOAA 9 451 
on 2001 May 9 with the Tenerife Infrared Polarimeter (tip; iMartmez Pillet et al.lll999l ) at- 
tached to the German Vacuum Tower Telescope in Tenerife, Canary Islands. The 0.5" x40" 
spectrograph slit was positioned across NOAA 9451 for approximately 70 min with no spa- 
tial scanning of the slit, w hile the solar image wa s kept stationary beneath the slit via a 
correlation tracking device (IBallesteros et al.l 1 19961 ). The mai n umbra from this dat aset was 
previously analyzed and presented as dataset 2 in the work of ICenteno et al.l (120061 ) . so only 
a brief description of the observational setup and format of the data is supplied here. Note 
the correct heliographic coordinates of the observed sunspot are S22°E20°. 

The spectral region obtained in these observations was recorded with a wavelength 
sampling of 31 mA pixel -1 and includes the photospheric Si 1 10827.09 A line, the upper- 
chromospheric He 1 10830 A multiplet (with triplet components at 10829.09 A, 10830.25 A, 
and 10830.34 A), and the telluric H 2 line at 10832.11 A. The tip instrument was used to 
record simultaneous spectral images of the four Stokes parameters (I, Q, U, V) for each of 
the 0.4" spatial pixels along the slit. Multiple images were coadded online to increase the 
signal-to-noise ratio, resulting in a final cadence of 2.1 s. 



After dark current subtraction, flat-field correction, polarization calibration, and re- 



- 3- 



moval of polarization cross talk (jColladoa 120031 ) . Stokes (/, Q, U, V ) profiles w e re inv erted 
separately for each line using the Milne-Eddington inversion code of lLagg et al.l (120041 ) . For 
the He I inversions an atmospheric model with one magnetic component was used, while a 
non-magnetic component was included in the Si I inversions to account for stray light. A 
continuum intensity space-time plot of the time series is given in Fig. [TJ alongside absolute 
magnetic field strengths, line-of-sight (LOS) velocities, and magnetic field inclinations in 
solar coordinates (see § 13.2.11 for more details) retrieved by the Si I and He I inversions. 

The Si I and He I LOS velocities retrieved from the inversions both show systematically 
decreasing velocity (i.e., blueshift) with time. This trend arises from the relative motion 
along the observer's LOS caused by the Earth's rotation, as these observations were obtained 
during local morning. Simple linear fits proved a suitable approximation to the trends and 
LOS velocity time series from each spatial pixel had these linear background trends removed 
prior to any form of temporal analysis. 



ANALYSIS METHOD 



Fourier phase difference analysis is a useful tool through which the propagation char- 
acteristics of waves may be determined. The form of an alysis used he r e is b ased on the 
standard Fourier equations that are discussed in depth by iKrijger et al.l (120011). This tech- 
nique has been used extensively i n the past for vario us solar studies (e.g., IJensen fc Orrall 
19631 . ISivaramanlll973l . lLiteslll984l . IWikstOl et al.ll2000l ) and it remains one of the most robust 
methodologies in use. 

The phase difference spectrum, A(f>(u), between two temporal signals measures the phase 
lag at discrete frequencies. For signals separated by some spatial distance, this lag is the 
cycle-time that it takes different frequency components to travel from the first location to 
that of the second. When waves propagate between two locations in a normally dispersive 
medium, phase difference spectra will show zero phase difference at low frequencies (where 
waves are evanescent and not propagating) followed by phase differences which increase in 
magnitude at higher frequencies (as high-frequency components travel more slowly than 
low-frequency components). 

The Fourier phase coherence spectrum, C 2 (z/), between two signals is a measure of the 
quality of phase difference variation. However, unless averaging in frequency is performed, 
the coherence between two Fourier components will be unity irrespective of the phase dif- 
ferences. The Fourier squared coherence of randomly distributed phase differences then 
approaches 1/n for averaging over n points in frequency. In this work coherence values are 
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calculated using an average over 5 frequency intervals, providing a "noise" level of 0.2 for 
randomly distributed (i.e., uncorrelated) phase differences. In the following phase difference 
diagrams, Fourier squared coherence is represented by the degree of symbol shading (white 
for 0; black for 1) while the symbol size represents the cross-spectral power - a measure of 
the co- variance between the Si I and He I LOS velocity signals. 



3.1. Trans-sunspot Wave 

If RP waves are due to waves propagating across the sunspot chromosphere, spectra 
calculated between the chromospheric LOS velocities from the umbra and those from pixels 
at sequentially greater distances into the penumbra should show phase difference values 
increasing linearly with frequency, becoming steeper with greater spatial separation of the 
signals. In this analysis, Fourier phase differences and coherences were calculated between 
the He I LOS velocity from pixel 55 (located in the umbra) and the He I LOS velocities from 
pixels at increasing distances into the penumbra. 

Figure [2] shows the output from such an analysis, where phase difference spectra be- 
tween the He I LOS velocity from pixel 55 and those pixel numbers listed in each panel are 
overplotted. For example, panel a contains A0 H e(t,55)^He(t,55) and A0 He (t,55)^He(t,54), while 
panel b contains A0 H c(t,55)^He(t,53) and A0 He (t,55)^He(t,52), where Ee(t,Y) denotes the tempo- 
ral He I LOS velocity signal from spatial pixel Y. Although groups of phase difference spectra 
are overplotted to increase the clarity of any relations, no clear form of propagating wave 
behaviour is seen. Note that around 4 mHz (i.e., 4-min period) values of increasing phase 
difference are observed when moving into the penumbra. However, the finding is marginal 
since the values of Fourier squared coherence in the spectra rapidly approach the "noise" 
level of 0.2 for randomized phase differences. 



3.2. Upward-propagating Waves 



Although a number of differing forms of propagating wave can exist in the outer atmo- 
sphere of a sunspot (e.g., fast/slow magneto-acoustic and Alfven waves) we shall restrict our 
analysis o f upward-propagating waves to that of field-aligned acoustic waves, as these were 
shown by ICenteno et al.l (120061 ) to describe the phase behaviour of 3-min waves in sunspot 
umbrae. The extension of these waves from travelling along near-vertical field lines in the 
umbra to travelling along inclined field lines which expand out over the penumbra is not 
unexpected since: 
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1. the magnetic field inclination increases smoothly from the umbral centre out through 
the penumbra (Figs. [T]d andCQg), 

2. the photospheric LOS velocity signals are fairly coherent across the umbra/penumbra 
boundary (Fig. [T]c). 

To study the possible propagation of field-aligned waves between two atmospheric 
heights we must first accurately determine the photospheric pixels that provide the lower 
atmospheric signal for the chromospheric pixels which lie above the penumbra. This is nec- 
essary because the field is significantly inclined here (Figs. [TJi and Wg) and velocity signals 
in the upper atmosphere will be spatially removed from their originating photospheric pix- 
els. The expected picture for field-aligned, upward-propagating waves is indicated in the 
schematic diagram of Fig. [31 where increasingly inclined field lines at the photosphere reach 
further into the chromospheric penumbra. 



3.2.1. Atmospheric Height Coupling 

In order to correctly pair spatial pixels between the photosphere and the chromosphere 
we require reliable determination of the magnetic field vector in solar coordinates. This 
is complicated by the 180° azimuthal ambiguity, whereby two equally valid but opposite 
azimuth orientations exist in the observer's coordinate frame. This uncertainty in the field 
azimuth impacts on the whole magnetic vector; the two differing azimuthal solutions yield 
different solar inclinations, 7'. 

To overcome this ambiguity we have implemented a "smoothest magnetic vector" form 
of ambiguity solution. Namely, a pixel region with the most realistic solution is selected 
as a trusted starting point (e.g., in the umbra where the true solution should be closest 
to vertical). Moving away from this seed region, either the 0° or 180° azimuth solution is 
chosen on a pixel-by-pixel basis to minimise the spatial variation in the three orthogonal 
components of the solar magnetic vector (B x , B y , B z ), where the ^-direction is normal to 
the solar surface. The field inclinations achieved in solar coordinates are shown in Figs. [TJi 
and for the inversion results from Si I and He I, respectively. 

Spatial pixels were then paired between the photosphere and chromosphere using the 
temporal averages of Si I inclinations in the solar coordinate framed, (7^)- Coupled with an 



1 Inclinations determined by the Sil inversion were chosen for this task because the greater signal-to-noise 
achieved in this line means that its magnetic vector, and hence solar inclination, is more reliably determined 
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expected height separation, AH, these inclinations provide pixel offsets between photospheric 
and chromospheric pixel pairs in the direction along the slit by AS = | tan (7^)] cos a AH/ s pix , 



where AH was taken as 1000 km (following the findings of ICenteno et al.l 120061 . for this 
sunspot umbra), s p ; x is the spatial sampling of the slit (^300 km pixel -1 ), and a is the angle 
between the field azimuth and the slit direction. Regions of pixels paired together in this 
work are outlined in panels b-g of Fig. [1] by dotted lines, while details of the pixel pairs 
are provided in Table [I] with their corresponding field inclinations. Although this approach 
uses the simplifying assumption that the magnetic field remains essentially linear between 
the two formation heights (i.e., there is no field curvature), it is somewhat justified by the 
resulting pixel pairs in Table [1] having inclinations that differ by ^7° - i.e., (7^) ~ (7h 6 ). 

The Fourier phase difference spectra resulting from these dual-height pixel pairs are 
presented in Fig. [H where spectra from groups of adjacent pixel pairs are again overplotted 
to enhance any relations. In contrast to the trans-sunspot case depicted in Fig. [2j the 
expected form of phase difference variation due to propagation (i.e., increasing values of 
phase difference with frequency) is clearly apparent in most of the panels. In addition, 
throughout panels a-g Fourier squared coherence values remain reasonably high. 



3.2.2. Dispersion Relation Comparison 



We make use of the equations provided in ICenteno et al.l (120061 ) which describe the dis- 
persion relation for vertical acoustic waves propagating in the presence of a vertical magnetic 
field within a stratified isothermal atmosphere with radiative cooling. The equations were 
modified for this work to simplistically mimic the first order effects that acoustic-like (low- 
(3 slow-mode) waves would experience when propagating along inclined field lines instead 
of purely vertically in a vertical magnetic field - i.e., cos 7' reduced gravity and I/COS7' 
increased path length. 

The solid curves in Fig. H] were calculated using the measured values of field inclination 
in solar coordinates along with the temperatu re (4000 K), radiative cooling time (55 s), and 
vertical height separation (1000 km) given by ICenteno et al.l (120061 ) for this sunspot umbra. 
Although only a simple approximation to the expected dispersion relation for such waves, 
these curves show an encouraging association with the measured data points. 



in comparison to that from the He 1 inversion. 
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DISCUSSION 



In this section we present our findings in the context of results from previous studies in 
an attempt to provide answers to a few of the outstanding issues surrounding the relationship 
that RP waves share with other forms of sunspot waves. 

It has been long k nown that sunsp ot chromospheres oscillate at differing periods in 
different spatial regions (jGiovanellil Il972l ) . Figure [5] displays the variation of Fourier power 
from the Si I and He I LOS velocities in the form of space-frequency diagrams. Individual 
power spectra from each spatial pixel have been normalized to the variance of the respective 
time series, resulting in white noise having power of 1 and 18.4 being the 99.99% signif- 
icance level of Poisson noise. Normalization was performed to aid in the comparison of 
spectral profiles between spatial regions that exhibit vastly different LOS velocity ampli- 
tudes. Ra ther than indicating some form of physical discontinuity (c.f., iTziotziou et al. 



20061 . 120071 ). the change from dominant chromospheric 3-min power to longer periods near 
the umbra/penumbra boundary in Fig. [5]& may just result from the magnetic field inclination 
becoming large enough to allow photospheric low- frequency (i.e., 5-min p-mode) power to 
tunnel through the high er-frequency acoustic cutoff (5.2 mHz) at the temperature minimum 
(IDe Pontieu et al.l I2004J ) . We note that the classical interpretation of an acoustic cutoff is 
effectively negated by the inclusion of radiative cooling in the modeled dispersion relation of 
§ 13.2.21 which allows wave reflection and transmission at all frequencies. However, the domi- 
nant chromospheric frequency in Fig. [5]o is modif ied by the magnetic field inclination, closely 
following the strong- field limit cos 7' relation of iBel fc Leroyi (119771 ). Power existing below 
the cutoff may be explained by the gradual transition from mainly reflected to transmitted 
waves around the cutoff (i.e., the slow turn-up in the curves of Fig. @J. Another possibilit y 
is the presence of unresolved structure in the chromospheric penumbra (jRiiedi et al.lll995l ). 
consisting of either two spatially-separated di stributions of fiel d inclination along the LOS 
or an uncombed magnetic field configuration (ILagg et al.l 120071 ). If the more vertical distri- 
bution has the measured field inclinations and the other has values ~20° larger, the acoustic 
cutoff curve in Fig. could be pulled to even lower frequencies. At larger field inclinations 
(i.e., further into the penumbra) little evidence is found of 3-min waves because power at 
5-min period vastly exceeds that at 3-min in the underlying photosphere (Fig. 0a) . 

The termination o f 3-mi n wav e patterns at the umb ra/penumbra boundary noted by 
Kobanov &: Makarchiki (120041 ) and iKobanov et al.l (120061 ) and the fact that not all 3-min 
wavefronts can be traced out from the umbra into the penumbra has been used to suggest 
that RP waves are not associated with similar waves in the umbra. A simple check for the 
linkage of either 3-min or 5-min waves between the chromospheric umbra and penumbra can 
be made by bandpass filtering the He I velocity time series. The spatial variation of He I LOS 
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velocities through the umbra and penumbra is presented in Fig. [6] before and after filtering 
in the period ranges 2.5 — 3.5 min and 4.5 — 5.5 min. It is clear that each of the 3-min umbral 
wavefronts has a rapidly diminishing counterpart in the penumbra (Fig. (H^), while each of 
the 5-min (i.e., RP) wavefronts has an only somewhat weaker counterpart within the umbra 
(Fig. [6]c). These findings once again support the picture of a continuous variation of RP 
wave behaviour through the entirety of the sunspot atmosphere. 

If RP waves are indeed the "visual pattern" of upward-propagating waves we expect that 
the wave velocity amplitude along the field would be essentially constant through the penum- 
bra, because each wavefront will have experienced the same degree of wave growth caused by 
the decrease of density with altitude. Examination of the unfiltered He I LOS velocity signal 
in Fig. \5)p, shows that the RMS LOS velocity decreases throughout the penumbra. However, 
the correct quantity to consider is the RMS velocity along the magnetic field vector. This 
was obtained from the RMS LOS velocity (solid curve in the right-hand panel of Fig. [6]a) 
using the measured field inclinations from the observer's LOS - resulting field-aligned RMS 
velocities are depicted in Fig. [6]a by a dotted curve. The field-aligned RMS velocities still 
show a decrease close to the umbra/penumbra boundary (from a greatly diminished contri- 
bution of transmitted power at 3-min period) but also a nearly constant value of ~1 km s~ x 
over pixels 45 to 25 of the penumbra, lending more credence to the "visual pattern" scenario. 



Our results support the conclusion of iKobanov et al.l (120061 ) that 3-min umbral waves 



are not the source of 5-min RP waves. However, we have additionally shown that they are 
in fact different manifestations of the same form of wave generated by a common source 
at the photosphere, their differences arising from the transmitted wave power available for 
propagation along differently inclined field lines. As such, the observed behaviour of waves 
in both umbrae and penumbrae can be explained without the need for abrupt changes in 
either density o r field orient ation at the umbra/penumbra boundary as postulated by, e.g., 



Tziotziou et al.l (120061 . 120071 ). 



5. CONCLUSIONS 

We have provided evidence that velocity signatures of RP waves observed in the He 1 10830 A 
multiplet are more compatible with upward-propagating waves than with trans-sunspot 
waves through careful consideration of the magnetic vector. Comparing the Fourier phase 
differences measured between paired pixels in the photosphere and chromosphere to the dis- 
persion relation for field-aligned acoustic waves, modified for inclined fields, points toward 
such waves (i.e., essentially low-/? slow modes) being responsible for the visual pattern. 
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Initially excited by a common source at the photosphere, waves experience increasing 
path length to the sampling height in the chromosphere with distance into the penumbra 
from travelling along increasingly in clin ed field lines - a scenario previously suggested by, e.g., 



Rouppe van der Voort et al.l (120031 ) and lBogdan fe Judge! (120061 ). For essentially constant (or 



weakly increasing) propagation velocities, delays of increasing magnitude will be observed 
in the arrival times of wavefronts at increasing radial distance through sunspot penumbrae. 
It is the pattern of delayed wavefronts that gives rise to the apparent outward motion of 
RP waves which may also explain the large range of observed wave speeds - the horizontal 
"speed" of the delayed wavefronts at the chromosphere depends on the rate at which the 
magnetic field inclines out through penumbrae, permitting either sub-sonic or super-sonic 
horizontal "speeds" for different magnetic geometries. This scenario also indicates that RP 
waves may occur at the edges of large pores since the existence of a penumbra is not necessary 
to support them; only sufficiently inclined field lines are required to direct the waves laterally. 

The German Vacuum Tower Telescope is operated on Tenerife by the Kiepenheuer 
Insitute in the Spanish Observatorio del Teide of the Instituto de Astrofisica de Canarias. 
The authors wish to extend their sincere thanks to R. Centeno, M. Collados, and J. Trujillo 
Bueno for providing this excellent data set for our analysis. 
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Fig. 1. — Space-time plots of: continuum intensity at 10825. 7±0. 3 A (a), absolute magnetic 
field strengths in Si I (6) and He I (e), LOS velocities in Si I (c) and He I (/), and magnetic 
inclinations in solar coordinates from Si I (d) and He I (g). Si I and He I velocities both 
have linear background trends removed; Si I velocities are scaled up by a factor of 8 to the 
dynamic range of the He I velocities. Upper solid white lines enclose a lightbridge in the 
umbra, while the lower white (black) line marks the umbral/penumbral (penumbral/quiet 
Sun) boundary. Temporal averages of the parameters (RMS values for LOS velocities) are 
shown in the right-most panels for both Si I and He I (solid and dotted curves, respectively). 
Regions of spatial pixels paired between Si I and He I (see § 13.2. ip are marked by horizontal 
dotted lines in panels b-g. Note oscillations in He I field strength and inclination are not real, 
but result from misfitting of the Stokes profiles associated with wave shocking. Observed 
shock profiles require two components; the single component inversion used here retrieves 
weakened, more inclined fields. However, the He I velocities retrieved still represent the 
general plasma motion. 
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Fig. 2. — Phase difference spectra between He I LOS velocity in the umbra and those at 
increasing distance into the penumbra. Symbol size and shading denote cross-spectral power 
and squared coherence, respectively. Panels show groups of phase difference spectra calcu- 
lated between umbral pixel 55 and the pixel numbers listed in the lower-left corners, moving 
from cases concerning pixels closest to the umbra/penumbra boundary (a) toward those in 
the middle penumbra (h). The Fourier squared coherence "noise" level has a value of 0.2 for 
randomly distributed phase differences in these data. 
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Fig. 3. — Cartoon schematic space-time diagram illustrating the form of pixel coupling be- 
tween the photosphere and chromosphere for the case of field-aligned, upward-propagating 
waves presented in § 13.21 Dark (light) grids denote umbral (penumbral) pixels, while dark 
(light) arrows indicate magnetic lines of force (i.e., wave paths) linking back to the photo- 
spheric umbra (penumbra). Note increasing delays in wavefront arrival time at the chromo- 
spheric sampling height because of increased propagation lengths along more inclined field 
lines. The horizontal and vertical axes are not to scale, resulting in magnetic field inclinations 
that appear different to the values actually retrieved from the Stokes inversions. 
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Fig. 4. — Same as Fig. El but for spatially-offset dual-height pairs of photospheric and 
chromospheric pixels. Curves show modified acoustic dispersion curves using the measured 
Si I field inclinations. 
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Fig. 5. — Space-frequency plots of Si I (a) and He I (b) velocity power. Spectra from each 
spatial pixel have been normalized to the variance of the corresponding time series; white 
noise has power of unity and 18.4 is the 99.99% significance level. Overlaid solid (dotted) 
curves are the acoustic cutoff modified by the Si I (Hel) inclinations, while horizontal lines 
mark the same boundaries as in Fig. [U 
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Fig. 6. — Space-time plots covering the sunspot umbra and penumbra, a) He I velocities 
retrieved by the inversion, b) He I velocities after bandpass filtering in the range 2.5—3.5 min. 
c) He I velocities after bandpass filtering in the range 4.5 — 5.5 min. Values in panels a and 
b are clipped to enhance wavefront visibility within the penumbra. White lines mark the 
umbra/penumbra boundary, while right-hand panels show undipped RMS LOS velocity. 
The dotted curve included in the upper RMS velocity panel shows the RMS velocity parallel 
to the field, after cosine correction for the inclination of the field from the observer's LOS. 
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Table 1: Spatial pairings between Si I and He I pixels and their corresponding magnetic field 
inclinations in solar coordinates 
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References. — 1 Panel of Fig. |4]in which the resulting Fourier phase difference spectra are overplotted 



